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Abstract 

A dual lattice vortex formulation of homogeneous turbulence is developed, within the Martin- 
Siggia-Rose field theoretical approach. It consists of a generalization of the usual dipole version 
of the Navier-Stokes equations, known to hold in the limit of vanishing external forcing. We 
investigate, as a straightforward application of our formalism, the dynamics of closed vortex tubes, 
randomly stirred at large length scales by gaussian stochastic forces. We find that besides the 
usual self-induced propagation, the vortex tube evolution may be effectively modeled through 
the introduction of an additional white-noise correlated velocity field background. The resulting 
phenomenological picture is closely related to observations previously reported from a wavelet 
decomposition analysis of turbulent flow configurations. 

PACS numbers: 47.27.Gs, 47.32. Cc 
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I. INTRODUCTION 



A considerable progress has been achieved along the last two decades concerning the 
kinematics of turbulent coherent structures, a fact intimately associated to the improving 
performance of computer and experimental resources 1]. However, the relevant dynamical 
properties of the evolution and interaction of the energy-containing eddies - believed by 
many to comprise the key for a fundamental understanding of intermittency and other 
turbulence characteristics - are still essentially unknown. As a concrete illustration of the 
present theoretical limitations, it is worth recalling the difficulties faced in the study of wall 
turbulence. Even though the main flow patterns have been identified in that situation 0,01) 
there is, for instance, no solid theoretical foundation for the logarithmic law of the wall. 

An ideal arena for the investigation of the dynamical and kinematical issues is provided 
by homogeneous isotropic turbulence. Direct numerical simulations have showed clearly 
that at moderately high Reynolds numbers the fiow is dominated by long-lived vortex tubes 
with small cross-sectional dimensions (defined around the Kolmogorov dissipation length) 
and sizes extending up to the integral scale 0, IE S 0, S]- also known, as set on a 

firm ground by Farge et al. 0], through wavelet decomposition methods, that most of the 
turbulent kinetic energy is carried by the vortex tubes, which are surrounded on their turn 
by a background incoherent fiow. 

Several analytical studies have addressed over the years the picture of homogenous tur- 
bulent fiows in terms of vortex tubes, either from the dynamical or kinematical viewpoints 
(see Ref. ^] for a comprehensive review). Among the former, a growing attention has 
been devoted to Lundgren's model based on the evolution of strained spiral vortices, 
which are transformed into tube-like structures and are probably generated in real fiows 
through shear layer instabilities [1,0]. In contrast, in the kinematical approach, the dy- 
namical details are bypassed, and an effective account of the statistical stationary regime of 
the "vortex tube gas" is attempted, as in the works of Chorin , where a connection with 
standard polymer statistical mechanics is drawn, and Hatakeyama and Kambe [16], whose 
focus relies on the properties of fiow configurations related to multifractal distributions of 
vortex filaments (modeled as Burgers vortices). 

Our initial aim in this paper is to establish, in Sec. II, an alternative formulation of 
the turbulence problem, incorporating into the usual stochastic approach the physical 
insight suggested from experimental and numerical investigations, which, as commented 
above, place vorticity coherent structures on a central stage. More specifically, we will 
implement, with the help of the Martin-Siggia-Rose functional formalism jl^, an exact 
statistical lattice vortex description of the fiow's dynamics jlEI3ll3' which contains, as a 
special case, the known dipole form of the Navier-Stokes equations. Next, in Sec. Ill, having 
in mind modeling matters (and, thus, non rigorous arguments), we use the lattice vortex 
formalism just developed to advance a phenomenological scheme describing the evolution of 
vortex tubes forced at large scales by stochastic forces. In particular, we also consider the 
effective force-force correlation function employed in the renormalization group analysis of 
turbulence [3, EO, 21) 2^1, decaying in Fourier space as k~^. In Sec. IV, we find that the 
stochastic perturbations due to the random external forcing may be effectively interpreted 
as resulting from the vortex tube advection by a white-noise correlated velocity background 
fiow. We determine the one- dimensional energy spectrum E{k) ~ k"^ of the background 
fiow, including its dependence upon the energy transfer rate, and the integral and viscous 
scales as well. It is interesting to note that a "thermal-like" energy spectrum, superimposed 
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to the Kolmogorov one, was indeed observed in the numerical wavelet analysis of turbulent 
configurations performed in Ref. To conclude, in Sec. V, we summarize and discuss our 
main results. 



II. DUAL LATTICE VORTEX FORMULATION 

As largely known, a systematic approach to the statistical description of homogenous 
isotropic turbulence, which is concerned with the flow's small scale properties, is yield by 
the stochastic generalization of the Navier-Stokes equations [17^ . 

dtVa + VpdpVa = yd'^Va - daP + fa , 

daVa = . (2.1) 

Above, fa = fa{x,t) denotes a gaussian random force, deflned at some large length scale L, 
with vanishing expectation value and the two-point correlation function 

{fa{x,t)fp{x',t')) = 6apF{\x-x'\)6{t~t') . (2.2) 



It follows from Novikov's theorem [2^ that energy is injected at large scales with pump- 
ing rate S = (foVn) = -^(0) = Dq. Furthermore, according to the standard Kolmogorov 
phenomenology ^], it is conjectured that dissipation takes place around the microscopic 
scale given by ~ Dq ^^^v^I^^ where viscous effects become relevant. The Reynolds number, 
depending only on the extreme scales L and r], is Re ~ {L/riY^^. 
Let the spatial part of force-force correlator be written as 

= Do exp(— m|x — x'l) , (2.3) 

with m = 1/L. An important feature of expression ()2.3j) . which does not necessarily hold 
for other admissible choices of the force-force correlation function, is that its inverse has a 
simple local form. Actually, we get, from (j2.3|) . 

F-Wx-x'\) = I {d^ -m^f6^{x-x') . (2.4) 

Notwithstanding the fact that the Fourier transform of F{\x — x'\) be regarded in principle 
as a regularized version of Dirac's delta-function in some appropriate functional space |25| . 
we will discuss, later on, dynamical effects related to the alternative definition 

F{\x-x'\) = — f d^kk-^ exp[ik ■ (f - x')] , (2.5) 



which has been a crucial ingredient in the renormalization group studies of turbulence [1 

23,|2l|,|23. In (EinD, the integration in Fourier space is bound to the region 1/L < k < l/rj. 



Observe that in this case, the mean energy input rate per octave is fixed to -Doln2, and we 
have F(0) = Doln{L/ri). 

Considering the bulk of experimental and numerical evidence that favors the picture of 
turbulence as a vortex tube gas, from now on our attention will be focused on the vorticity 



3 



dynamics implied by the stochastic Navier-Stokes equations. It is convenient, thus, to work 
with the stochastic Helmholtz equation, straightforwardly derived from Eqs. ()2.1|) as 



dtUJa + Vfsdf^Ua - UJpOpVo, = 1^8 UJ^. + f* , (2.6) 

where /* = eap-ydpf-^ and = ^ap-tdpV'^ is the vorticity field. We obtain, using (|2.2p . the 
correlator 

{ax,t)fi{x',t'))^D^p{\x-x'm-t') 

= eapa(^l3'yv9pd!y{faix,t)fr,{x',t')) 

= {d^dp-5apd^)F{\x-x'\)5it-t') . (2.7) 

We state now, in field theoretical language, what is meant by the stochastic evolution 
problem. Defining uj^{x) as the vorticity field at a certain time instant Iq, we are interested 
to find the probability density functional Z = Z[uJa{x) ,ti\u^{x) ,tQ] for the observation of 
vorticity uja{x) at a latter time instant ti. Within the path-integral version of the Martin- 
Siggia-Rose formalism [3], it follows, from Eqs. (j2.6|) and (j2.7|) . that 

Z =U j DuJaDuJaexp{iS) , (2.8) 



where A/" is a normalization constant [26[, to assure that / Dua{x)Z[uJa{x),ti\u^{x),tQ] = 1, 
and 



S = [ dt{ I d^xCja{dtUJa + Vj^dj^UJa — UJpdpVa — vd'^UJa) 
Jto J 

+ i [ d^xd^x'uJaix,t)Daf3i\x- x'\)LUf3{x',t)} . (2.9) 



An effective model of turbulent dynamics would be naturally attained if the velocity and 
vorticity fields that appear in ()2.9|) were expressed as a sum over the contributions produced 
exclusively by relevant fiow profiles. The basic difficulty here regards the selection and 
parametrization of such configurations. A promising starting point is to establish a set of 
"building blocks" that could be used to represent the usually observed coherent structures. 
We recall that vortex sheets or tubes, in particular, can be exactly obtained in a simple way as 
linear combinations of elementary closed vorticity rings, through a lattice vortex construction 
"3 3 originally devised in the realm of superfiuid physics jl^]. Just define a cubic lattice. 



with spacing parameter e ^ (i.e., much smaller than the Kolmogorov dissipation length 
rj), whose sites are written as Xp = e{piXi + P2X2 +^32:3), where the pi's are integers. The 
vector position Xp is taken to be the common vertex of three plaquettes oriented according 
to the unit vectors x^, as shown in Fig. 1. In a self-evident notation, an arbitrary plaquette 
is completely charaterized by the vector doublet V = {xp,Xcr). Furthermore, by definition, 
the plaquette's boundary dV is identified to a line vortex (vortex tube with vanishing cross 
section) which carries vorticity fiux 4>a{xp, t). Of course, a square line vortex is an ill-defined 
mathematical object, due to the divergence of the velocity field on its corners 113] . However, 
as a simple regularization procedure, we impose an ultraviolet cutoff A = 1/e in the Fourier- 
transformed kernel of the operator that maps vorticity into velocity. In rough terms, this 
is equivalent to replace the line vortex by a vortex tube with a cross section of radius ~ e. 
It is important to remark that the lattice vortex is in fact an "overcomplete" basis for the 
description of general flow conflgurations. 
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Considering the plaquette V = {xp,Xa), let Xa{s) G dV be a point parametrized by the 
arclength < s < 4e of the oriented boundary path that starts at the reference point Xp 
and ends at Xa{s). The vorticity field associated with this plaquette is 



(x;q = 0<^(fp,t)5(ni)5(n2)^a;Q(s) , (2.10) 



where rii and n2 are the coordinates along the normal and binormal directions on the line 
vortex (the binormal vector is defined as Xa)- The central idea underlying the lattice vortex 
representation is, then, to substitute ()2.10|) in the Martin-Siggia-Rose action ()2.9|) and per- 
form afterwards the sum over all the plaquettes. Introducing = daCjj^ — dpuia, we get, 
for a single plaquette, the following relations: 

J Jav Jav 

d^XUa{VpdpUJa - Upd/sVa) = <j)a{Xp, t) f dx^F^fiVp , 

Jav 

d^XUJad'^LOa = 'PcriXpjt) f dXaO'^UJa = —(j)a(Xp,t) f dXadaFaB . (2-11) 

Jar Jav 
Thus, the Martin-Siggia-Rose action becomes 

S = f dtiY^f dXa[dt(f)ad~'^dpFo,p + (j)„Fo,pVp + u(f)^dpFo,p] 
Jto p Jap 

+ d'xd'x'F^f,{x,t)F{\x- x'\)F^f,{x',t)} . (2.12) 

It is worth noting that the action ()2.12|) is invariant under the local transformation — ^ 
il^a + daX- III fact, it turns out that the transition probability for small time intervals, derived 
from ()2.8|) . with ()2.12p . is well approximated by the expectation value of a product of loop 
operators, computed in a non-local, three-dimensional, U{1) gauge theory. 

Drawning upon the gauge field theory correspondence, we define now the dual field 
strength, 

(l>a{x,t) = -eap^Fp^{x,t) , (2.13) 
which satisfies to da(t)a = and, additionally, 

Fal3 = ; 
-FapFafi = 4>a4>a , 

dXaF^pVp = eaf3ydXaVf3^^ . (2.14) 

We find, substituting the above relations in ()2.12p . 

S = - [ dt{J2 ^ dXaeaPy[dt(f)a{Xp,t)d~'^dp(j)^ + (j)„{Xp,t)(f)pV^ 

Jto -p Jav 

- u(l)^{xp,t)di3(f)y]+i I d^xd^x'(f)a{x,t)F{\x - x'\)(f)a{x',t)} . (2.15) 
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The solenoidal constraint for the dual field can be imposed in the path-integration ()2.8|) 
by means of an auxiliary scalar field A, which is nothing but a Lagrange multiplier. More 
concretely, we take 

Z = U j D^aD(l)^D\exp{iS) , (2.16) 



with 

S = j dt{J d^x{Laix,t) + daX)(j)f3ix,t) 

+ i J d^xd^x4)a{x,t)F{\x - x\)^a{x',t)} . (2.17) 
In (|2.17p . the whole dependence on the 0-fields is implicit in the non-local "source" term 
L„(x,t) = -Y.f dx' fa,i3^[dt(pa{xp,t)d~^dp 

P JdV 

+ (p„{xp,t)vp{x' ,t) - u(p„{xp,t)dp]5^{x' - x) . (2.18) 
To proceed, we define the Fourier transform of La{x,t), 

La{k,t) = / d'^xexp{—ik ■ x)La{x,t) = iy^^ f dx'^eajs-y exp{—ik ■ x') 
x[^dt(f)a{xp,t) + i(f)^{xp,t)vp{x',t) + u(j)^{xp,t)kf:i] . (2.19) 

— * 

Writing, for a given plaquette V = {xp,Xa), the boundary position vector as = Xp + ^, 
with Va{x') ~ Va{xp) + ^r]dr^Va{xp) , wc obtain 

La{k,t) = i'^expi^-ik ■ Xp)eaf3-y{gay{k)[^dt(l)a{Xp,t) 

p 1^ 
+z0^(fp, t)vp{xp, t) + vkp(j)„{xp, t)] + 5'a7,»y(^)0a(^p, t)d.^vp{xp, t)} , (2.20) 



with 

da-yik) = f d^yexp{-ik ■ i) , 

J dV 

9a-y,vi^) = -i £^ d^^iri exp{-ik ■ = '^9a-y{k) . (2.21) 

In the limit e ^ 0, keeping k <^ 1/e, we have, asymptotically, gap{k) = ie'^ea/S'yk^, and, 
therefore, 

La{k,t) = e^^exp(-zA; ■ Xp){nQ^(A;)[9t0^(fp,t) + uk^(f)f3{xp,t)] 
p 

+i(f)p{Xp,t)[5af3V^{Xp,t)k^ - idaVf3{Xp,t) - kaVi3iXp,t)]} , (2.22) 

where we used Uaf^^k) = 5ap — kakf^/k"^, the Fourier-transformed projector on transverse 
modes. The continuum limit of the above sum is defined through the substitutions 

^ X , ^ ^ 1 d^f , 0^ ^ . (2.23) 



We find 

L^{k,t) = fl^p{k)[dtMk,t) + ueMk,t)] + na{k,t) , (2.24) 

where 



(j)a{k,t) = J d^x ex.Y>{—ik ■ x)(j)a{x,t) , 

Vlaikit) = / d^x exp{—ik ■ x)v(}{di3(f)a — da4>f3) ■ (2.25) 



There is a simple connection between 0q, and the velocity field Va- Taking the Fourier 
transform of the vorticity field, we find 

i,<.(fc,() = /d'fexp(-«.fK(f) 
= d.„exp(-«.fM,(x„i) 

^ JdV 

= '^gpa{k)exp{-ik ■ Xp)(f)i3{xp,t) 
p 

= ie^eap^k/3^exp{-ik ■ Xp)(j)^{xp,t) . (2.26) 
p 

Recalling ()2.23p . we get Ua = ^ai3-tdi3<t''y i^i the continuum limit. Since = ^a/S'ydpv^, we 
immediately conclude that the fields Va and (pa differ only by a gradient, which means that 
Va = Ilai3(f>i3- As a cousequeuce, if ()2.22|) is taken back to real space, we get 

La{x, t) = dtVa - yd'^Va + Vp{dp(j)a - 9„0/3) . (2.27) 

We define at this point the additional scalar field C = [dt—vd'^]~^X, and impose (pa = Va+daC, 
so that 

La + daX = dt(pa - vd'^(pa + Vf3{dp(pa - OaCpf^) ■ (2.28) 

Since the action ()2.17|) is quadratic in (pa, it is possible to evaluate the exact path- integration 
over the dual fields. Using ()2.16|) . ()2.17|) and ()2.28p . the result is an effective (and exact) 
expression for the probability density functional Z, 

Z =M j D(paD\ exp(z5^) , (2.29) 

where 

= -^j^ dt J d^xd^x[dt(pa- l'd'^(pa+ Vp{dp(pa- da(pf3)\x 

X F-\X - x')[dt(Pa - l^d^<Pa + Vp{dp(Pa ' da<Pp)].' ■ (2.30) 

The field theory given by ()2.30|) may be obtained directly, along the Martin-Siggia-Rose 
formalism, from the stochastic differential equation 

dtCpa + Vp{dp(pa - da(pp) = l'd'^(pa + /a • (2.31) 

For vanishing external forc es, t he above expression reduces to the usual dipole version of 
the Navier-Stokes equation [l5|. Substituting (pa by Va + daC in (j2.31|) . we get the original 
stochastic Navier-Stokes Eqs. 1)2.11) . with pressure given by P = dtC — vd'^C, — . 
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Does Eq. ()2.31|) yield any advantage over the standard Navier-Stokes formulation? Di- 
rect numerical simulations based on ()2.3H) would probably have the same computational 
cost than the ones which usually rely on the Navier-Stokes equations, since both versions 
involve at least two Fourier transformations per iteration cycle. In practice, the above de- 



scription provides an alternative approach to large eddy simulations [IJ], or the analysis of 



phenomenological aspects of vortex tube dynamics, as put forward in the following consid- 
erations. 



III. STOCHASTIC VORTEX TUBE EVOLUTION 

We are now interested to investigate the evolution of a closed vortex tube F, with small 
linear cross-sectional dimensions (of the order of rj) and subject to the action of large scale 
gaussian random forces. In a first approximation, we regard the tube as a vorticity filament, 
parametrized by the curve Xa = Xa{s,t), and carrying total vorticity flux 0. The vorticity 
field is given by 

Ua = (pS{ni)6{n2)—Xa{s,t) , (3.1) 
as 

where, similarly to the former plaquette's definitions, ni and n2 indicate the normal and 
binormal coordinates along the line vortex. 

The assumptions taken in ()3.1|) that the vorticity flux is time-independent and that 
cross-section fluctuations may be neglected are imposed as phenomenological constraints. 
Our results will be expected to hold to the extent that phenomena like vortex breakdown, 
vortex merging, etc. do not affect the vortex tube evolution. Such flow regimes have been 
well verified in the numerical and real experiments where vortices are mostly advected by 
the background flow, during their mean life-time, in agreement with the flux conservation 
Kelvin theorem. This state of affairs gives in fact the physical basis that supports the 
somewhat popular choice of modeling vortex tubes by means of Burgers vortices, or similar 
configurations. 

Our first task here is to apply the information provided by (j3.H) in the effective action 
(j2.3(J|) . In the limit of vanishing viscosity, we are left, therefore, with the evaluation of dt(j)a 
and Vp{df3(f)a — da<t>i3). The latter quantity is just minus the Lamb vector. In fact, using 

= ^a/3'rd/3(p'y a straigthforward computation leads to 

Vp{di3(f)a - da(t)f}) = eafi-yUJpV^ ■ (3.2) 

On the other hand, to find dt(f)a, let us imagine as an auxiliary construction that the line 
vortex is advected by a divergence- free field ^q(x, t), defined on all space, and which satisfies 
the boundary condition ^a(x, t) = Xa(s,t) on F. We have, then, |28| 

dtu = V X {^x u) . (3.3) 

Observing that dtUJa = <9t[eQ,/3^9^0^] = eap-ydpOtcj)^, we get, from (j3.3p . 

dt<Pa = ea^yXpu^ + daX , (3.4) 

where A is an arbitrary field. We are ready to substitute ()3.2|) and ()3.4|) into ()2.30|) . Intro- 
ducing 

= eo,i3^i)p{s,t)—x^{s,t) , (3.5) 
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where ifja = Xa — Va, we obtain 

•S*,^ = S^^ + Sxii) + Sxx , (3-6) 



with 



S^^ = ^ dt ds ds'iP^{s,t)F-\x{s)-x{s'))^Pi{s\t) 



4 Jto 



2 Jto 
i 

4: Jto 



Sxx = 7 r dt I d^x I d^xdaX{x,t)F-\x - x')daX{x,t) , (3.7) 



where p{t) is the length of the vorticity filament. The integration over A gives 

Z = Af I D^^ exp{iS^) , (3.8) 



S^ = -^ dt ds ds'ij^{s,t)Ua,pF~\x{s) -x{s'))ij^{s',t) . (3.9) 



where 



4 Jto 

Note that while 5*^ is a functional of tp^, the projection of tpa on dxa/ds (that is, the longi- 
tudinal component of ipa) maps the line vortex into itself. The singularities that eventually 
appear in integrand of ()3.9|) may be circumvented in a physical way, replacing the original 
vortex filament by a vortex tube, through the substitutions 

ds — > d^x , ds' — > d'^iJ" , 

F-^(f (s) - f(s')) ^ F-^(f - x') , (3.10) 
where x = {rii, n2,s), x' = {n[, n'g, s'), and 

1 1 

h{ni, n2) = — 5- exp[ — -{nl + nl)] . (3-11) 

We assume that the curvature radius of the vortex tube is much larger than the Kolmogorov 
dissipation length, an hypothesis supported by observations. A necessary condition for this 
property to be preserved in time is 

V\ds^^\ < \^^\ . (3.12) 



In practical computations, we may work with a straight vortex tube, taking s = z, ni = x, 
and n2 = y (one can figure it out as a circular vortex tube with infinite curvature radius). 
We deal below with two specific examples of external stochastic forcing, given by ()2.3j) and 
fl2.5|) . which will be named models A and B, respectively. A more concise expression for S^, 
compared to ()3.9|) . follows in general, relying basically on the slender profile of the vortex 
tube. 
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Analysis of Model A. 

To evaluate it is necessary to consider the operator kernel, 



U^pF-\x -x') = ^-f^{6^/3 - d-^d^dp){d^ - m'fd^x - x') . (3.13) 

If this expression is substituted into ()3.9|) . considering ()3.10|1 and ()3.11|) . a number of terms 
is obtained, hierarchically organized according to the powers of the dissipation length 77 — >■ 
defined in their coeficients. We will retain in the expression for only the dominant term, 
corresponding to the smallest power of 77. Using rotation invariance around the z axis, and 
neglecting derivatives of ip^ along the z direction, as it follows from ()3.12|1 . this prescription 
effectively amounts to perform in ()3.9|) the replacement 

Yi^pF-\x-x') ^ {dlf5\x-5') , (3.14) 

where {d±Y = d1 + dl. We obtain 



= ^TTTT^ / / ds[i,i{s,t)f . (3.15) 



Although (|3.15|) is an apparently simple quadratic action, the time-dependent spatial inte- 
gration limit p{t) renders the analytical evaluation of Z difficult. Nevertheless, the problem 
looks amenable of numerical investigation through the use of Langevin techniques j22| . 

Since we are discussing the time evolution of a vortex tube, the relevant physical question 
one may ask is concerned with the probability density functional of finding the tube in a 
certain geometrical configuration. In a first instance, this seems to be an intricate problem, 
once any individual vortex tube "world-line" to be considered in the path integration is 
accounted for by a large number of configurations ofipa- A simple solution of this degeneracy 
problem may be obtained, however, by means of the "minimal mapping" depicted in 
Fig. 2. The essential idea is to keep track of the vortex tube evolution for a very small time 
interval 6. We decompose the time evolution in two steps. First, the tube T{t) is mapped 
into r*{t) through its self-induced velocity field Va- Next, the stochastic perturbation 1/^^ 
takes T*(t) to the final configuration T(t + 6). The mapping sequence is Xa ^ x'^ ~^ ^'L 
with 

Xq, Xq -|- 5f q, , 

x'^ + 5iIjc,. (3.16) 



Let 7 be the plane that contains x'^ and is normal to r*(t). Then, is just the vector 
parallel to 7 that connects x'^ to the vortex tube T{t + 5). We have 

i^i{s,t) = i^l{s + 5Us.t),t) + 0{5^) 

= <(s,t) + 5^,9,<(s,t) + 0(52) , (3.17) 

where ips = ipadxa/ds. The expansion (|3.17|) implies that ip-^ ~ "^a = 0(5), and so may 
be substituted by in (j3.15|) . We find, thus, that the probability density functional for 
the transition r(to) r(ti) of the vortex tube configuration may be defined as 
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Clearly, the original field degeneracy is removed, and there is in ()3.18p an one-to-one corre- 
spondence between the vortex tube integration paths and the fields '?/'°(s,t). 

Analysis of Model B. 

The computational steps are exactly the same as the ones performed in the former case. 
The only technical difference is that the analogous of Eq. (j3.13p is written now in Fourier 
space: 

IlapF-\k) = (S^^e - Kkpk) . (3.19) 

The dominant contribution to ()3.9p . of order l/?7^, comes from the substitution 

Iio.pF~\k) ^ • (3.20) 

We get, similarly to (|3.18|) . the probability density functional 

Zs=AfJ D^IjI exp{-^-^ dt ds[^l{s, t)f} . (3.21) 

A remarkable feature of model B, as it may be easily inferred from ()3.21|1 . is that there is 
no dependence of the probability density functional Zb upon the integral scale L = 1/m (as 
it occurs in model A, for instance). 

In order to establish a connection between the above models and observed features of 
turbulent fiows, a slight modification of expressions ()3.18|) and 1)3.211) is necessary. In prin- 
ciple, the Martin-Siggia-Rose framework implemented by ()2.29|) and ()2.30|) is expected to 
provide a bona fide statistical modeling of vortex tube motion if an ultraviolet cutoff appears 
dynamically at a frequency |a;| ~ 1/^?,, where t.^ ~ r]'^^^ is the eddy turnover time at the 
Kolmogorov length scale. The simplest way to find improved versions of ()3.18j) and (|3.21|) . 
thus, is to replace the Dirac's delta factor in ()2.2|1 by a regularized expression like 

6R{t-t') = ^expi-t-'\t-t'\) , (3.22) 

and relax the cutoff prescription for the field ipa in frequency space. As a consequence, 
if all the steps leading to ()3.18|) and ()3.21|) are evaluated again, taking into account the 
modifications due to ()3.22j) . we will get, for both models A and B, the general result 

D<exp{-c y dt ds[{t^dt^lf + «)2]} , (3.23) 

where the constant c = c(m, r], Dq) is defined in the same way as before. The spatial 
wavenumber cutoff is hidden in ()3.23p . insofar as it will not have any relevant role in the 
forthcoming arguments [the cutoff is much smaller than k^j ~ l/r], according to ()3.12p ]. 

IV. BACKGROUND VELOCITY FLUCTUATIONS 

It is interesting to note that the probability density functional ()3.23|) is completely equiv- 
alent to the one derived for the problem of random advection of a vortex tube by a back- 
ground velocity field. In this way, we can draw a correspondence between the former effective 
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description, based on the analysis of the stochastic Navier-Stokes equations, and reahstic 
properties of turbulent flows. In more precise words, let Va{x,t) be the velocity of the back- 
ground flow, which is assumed to be a random gaussian fluctuating field, with vanishing 
mean value and correlator 

Mx, t)vp{x\ t')) = gIlo.p5\x - x')6R{t - t') . (4.1) 

It follows, therefore, that the one- dimensional background energy spectrum is given by 

and that the path-integral expression (|3.23p holds, with 

c=f!. (4.3) 

It is straightforward to prove ()4.2|) from the Fourier transform of the velocity-velocity corre- 
lator (j4.ip . Below, we discuss in more detail how (j3.23|) arises from (j4.H) . with the specific 
parameter definition (j4.3p . Other velocity correlators would work as well. However, (j4.H) is 
particularly attractive in view of its direct relation to numerical observations 9]. 

The probability density functional to have a certain background velocity field Va{x^t) in 
the region VLt enclosed by a vortex tube, for the time interval to < t < ti, may be written as 

V = {UijSiv^ixi, tj) - Va{xi, tj))) , (4.4) 

where {xi,tj) denotes a discretized space-time position defined in the set of world-hnes gen- 
erated by the vortex tube evolution. Using the Fourier representation of the delta function, 
Eq. fl4.4|) becomes, in the continuum limit, 

V = J\f f D^o.exp{i f'dtf d^x^aVa){exp{-i f ^ dt [ d^x^^^Va)) ■ (4.5) 

J Jto Jflt Jto Jvtt 

Resorting to the gaussian random behavior of the background velocity field, we are able to 
compute the above expectation value. Using 1)4.11) we find, 

V = U f D^^exp{i r dt j c/^f^^-y^) exp[-| f^dt dt' 

X / d^xUx,t)6R{t-t')U^p^p{x,t')] . (4.6) 
J fti 

Since Va = HaisV/s, we may integrate over the field to get 

V cx exp{-;^ dt I d'xlitr^dtv^f + [v^Y]} . (4.7) 

If the vortex tube has a small circular cross section of area rrrj'^, we can replace /^^ d^x by 

Tvrf Jq^^^ ds in ()4.7|) . Furthermore, to find the transition probability density functional Z 
for the vortex tube evolution between configurations Tito) and r(ti), we (i) decompose the 
velocity field in transverse and longitudinal components to the vortex tube tangent vector, 
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viz., Va = v-^ + v^ai (ii) integrate over the longitudinal components w^, and (iii) introduce the 
"minimal velocity field" f° in close analogy with the previous definition of if)^. We obtain 

Z=U j Dv'iexpi-^ l"^'^ ds[{t,d,vy + ivy]} . (4.8) 

Therefore, identifying f ° to ip^, we have just found ()3.23p again, with c given by ()4.3|) . From 
fl4.2|) and ()4.3|) . we can predict the form of the one- dimensional energy spectrum for models 
A and B (disregarding numerical prefactors): 

E.(^) - . (4.9) 

It is useful to compare the Kolmogorov's spectrum Exik) ~ D^^^k^^^^ with the above 
expressions. We may estimate, relying on Kolmogorov phenomenology, that (p ~ DH^rj^^^ 
and ~ D^^^^rf^^. At the dissipative wavenumber A;^ ~ l/?7, we define the Reynolds 
number dependent dimensionless ratio 

(4.10) 

where E(kr^) is the background spectrum for a given model. It turns that for model A, we 
get a = —1 while for model B, a = 0. More generally, it is not difficult to realize that the 
family of gaussian stochastic forces described by 

Fik) ^ + m^)-^ , (4.11) 

with /? > 3/2 leads to (pTTTHl with a = 3-2/3. 

A numerical wavelet analysis by Farge et al. of the direct numerical simulations 
carried out by Vincent and Meneguzzi at moderately high Reynolds numbers, reveals the 
existence of a background fc^ one-dimensional energy spectrum. The turbulent fiow may be 
depicted as a vortex tube gas surrounded by incoherent fiuctuations, the latter having their 
kinetic energy equiprobably distributed over the spatial Fourier modes. It has been suggested 
in Ref. |^ that the dissipation at the bottom of the inertial range would be preceded at 
larger scales by some coherent-to-incoherent energy transfer from the vortex tubes to the 
background field. A fraction of the vortex tubes would be disrupted in a conservative way, 
so that the transformation of their mechanical energy into heat would occur afterwards from 
the background fiow. One may conjecture that the integral length scale is irrelevant in this 
sequence of small scale events. In that case, we have a = 0, as in model B, which is actually 
the scenario indicated by the numerical results, where Q — 0.1 for the Taylor-scale Reynolds 
number Rx = 150 (equivalent to Re — 10"^, according to Lohse [s^). 

V. CONCLUSION 

We investigated in this work both formal and phenomenological aspects of homogenous 
isotropic turbulence, within the stochastic modeling of vorticity dynamics. A rigorous sta- 
tistical lattice vortex description of turbulent fiows was established, yielding the basis for a 
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subsequent phenomenological discussion of the problem of the random evolution of vortex 
tubes, commonly observed in experiments and numerical simulations. Since the advection 
of vorticity coherent structures is ultimately caused by the background flow, according to 
Kelvin's theorem, we interpret the stochastic method as an effective tool for computing 
the evolution of vortex tubes. We were able to find in this way a plausible form for the 
background velocity-velocity correlator, and, as an immediate consequence, the background 
one-dimensional energy spectrum. We found a satisfactory agreement with the recent nu- 
merical analysis of Farge et al. 9], where a thermal-like spectrum was clearly noticed for 
the background flow. In particular, we observed that the gaussian correlator ()2.5p . used in 
the renormalization group approach to turbulence 3, 2^, 1^ 2^1 is likely the correct choice 



(model B of Sec. IV) for the derivation of phenomenologically meaningful results. Further- 
more, it would be important to improve the connection between the stochastic modeling 
and the numerical results concerned with anisotropic effects, as the reported zero helicity 
distribution peak for the incoherent fluctuations 

There is a strong numerical evidence that the vortex tube gas accounts on its own for the 
Kolmogorov's spectrum 31, 3^. Regarding the background flow, our analysis suggests 

that it has a twofold character, involving the combination of the "eddy noise" 2^1 forcing, 
effectively modeled by ()2.5|) . and of configurations which satisfy the energy equipartion 
principle. The picture that emerges - to be explored in further analytical and numerical 
works - is that these two facets of the background fluctuations are self-consistently related 
to the vorticity coherent structures. While the force- force correlator ()4.11|) . with f3 > 3/2 
is a reasonable choice for a rigorous study of the turbulence problem, it becomes useless 
when considered in the simplified phenomenological perspective addressed in Sec. IV. On 
the other hand, model B is favored by the force of numerical observations, since it copes well 
with the tripartite phenomenological stage set up by the vortex tube gas, stochastic eddy 
noise, and the thermal-like background flow. 
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FIGURE CAPTIONS 



FIG. 1: The three oriented plaquettes which have the common reference position Xp, and carry, 
on their boundaries, vorticity fluxes $i, $2, and $3. 

FIG. 2: The vortex tube r(i) evolves, during the small time interval 6, to the new configuration 
T{t + S). The intermediate dashed tube T*{t) corresponds to the transport provided by the self- 
induced velocity field Va- 
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